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ABSTRACT 


This report presents the techniques applied to MSS- 
LANDSAT data in the land-use mapping of Criciuma region (Santa Catari 
na state , Brazil) and the results of a classification accuracy 
estimate tested on the resulting map. The MSS-LANDSAT data digital 
processing involves noise suppression, feature selection and a hybrid 
classifier. The accuracy test is made through comparisons with aerial 
photographies of sampled points. The utilization of digital processing 
to map the classes agricultural lands, forest lands and urban areas is 
recommended, while the coal refuse areas should be mapped visually. 
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1. INTRODUCE 


0 objetivo deste relatorio e apresentar os resultados da 
avaliagao do desempenho das tecnicas de analise digital de dados MSS- 
LANDSAT aplicadas ao mapeamento do uso do solo na regiao carbonTfera de 
Santa Catarina. 

A area de estudo para esta avaliagao corresponde ao modu 

k 

lo de Criciuma, previamente definido em Valeriano e Pereira (1983). 

Todos os procedi men tos de tratamento digital de imagens 
foram realizados no analisador de imagens IMAGE-100 do INPE. Os algorit 
mos utilizados estao descritos em Ribeiro et alii (s.d.). 

Apresenta-se a metodologia de processamento digital de 
imagens MSS-LANDSAT aplicada a area de estudo, bem comoo mapa resultani 
te deste procedimento. Descrevem-se as tecnicas de interpretacao visual 
de fotografias alreas utilizadas na elaboragao do mapa de verdade ter 
restre, o qua! tambem e apresentado neste relatorio. 

Os criterios e a metodologia para a estimativa da exatj_ 
dao de cl assifi cacao do mapa de uso da terra derivado de dados MSS- 
LANDSAT encontram-se neste relatorio. 

2. CLASSIFICAQAO DIGITAL DE DADOS MSS-LANDSAT 

2.1 - prE-processamento 


Ao realizar urn tratamento digital de imagens com fins de 
classificacao automatica, deve-se, prel iminarmente, eliminar as informa 
Coes espurias da cena como ruidos da aparelhagem e interference da 
atmosfera. 


Para eliminar ruidos causados por falhas de aquisicao ou 
gravacao dos dados, aplicou-se no canal MSS 6 (contendo uma faixa em 
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branco) o algoritmo DRUIDO que, ao localizar tal faixaderuTdo , a subs, 
titui pel a media das faixas adjacentes. 

0 efeito de "stripping" decorrente da sensibilidade dif£ 
rencial dos seis sensores de cada banda espectral de aquisigao de da. 
dos do MSS-LANDSAT foi atenuado pelo algoritmo CORRECAO RADIOMETRICA, 
desenvolvido para este fim. 

* 

A eliminagao da radiancia retroespalhada pel a atmosfera 
foi realizada pel a subtragao, em toda cena, dos valores medios de cada 
canal do MSS-LANDSAT, obtidos de "pixels" local izados em areas de som 
bra de relevo, 

2.2 - GERACAO DE NOVOS ATRIBUTOS 


Dentro da terminologia estatTstica utilizada em classify 
cacao, atrihuto e uma caractenstica mensuravel do conjunto que se pre. 
tende classificar. As radiancias espectrais do "pixel" imageaao, tradj£ 
zidas em valores de LANDSAT (LANDSAT "number") em cada canal, sao os 
atributos originais oferecidos pelo sistema MSS-LANDSAT. Entretanto, 0 £ 
tros atributos podem ser gerados a parti r dos originais para a classj. 
fi cacao da area. 

As razoes entre canais fornecem rapidamente (dada a facj. 
lidade de operagao) informagao sobre a diferenga das radiancias espec 
trais do alvo naqueles dois canais. Alem disto, atenuam o efeito de ga 
nho causado pela orientagao do terreno em relagao ao Sol (efeito de 
iluminagao e sombra). 

Dos quatro canais MSS-LANDSAT foram obtidas imagens do 
resultado da divisao de cada possTvel par de combinagao. Estas imagens 
foram denominadas de acordo com as seguintes razoes: 


RAT 5/4. 
RAT 6/4. 
RAT 7/4. 


RAT 6/5. 
RAT 7/5. 
RAT 7/6. 
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2.3 - SELEQSO DE ATRIBUTOS 


Nos procedimentos de cl assifi cacao automatica, o analis^ 
dor de imagens IMAGE-100 comporta apenas 4 canais. Portanto, dos 10 atri_ 
butos disponTveis (6 razoes entre canais e os 4 originais) foram sel£ 
cionados os 4 que melhor diferenciam os aspectos do terreno a ser cla_s 
sificados. 

* 

0 algoritmo SELECAO DE ATRIBUTOS calcula a distancia J.M.* 
entre classes caracteirizadas por treinamento em cada combinagao de qua 
tro atributos. Em seguida, seleciona a melhor destas combinagoes atra. 
ves dos criterios de maior distancia minima e maior distancia media eji 
contrada. 


A probabil idade de erro de classificagao, baseada nas di£ 
tribuigoes fornecidas pelas amostras de treinamento s pode ser inferida 
dos limites maximo e mTnimo oriundos da distancia J.M. 

As classes de uso da terra de maior interesse para este 
trabalho (rejeito piritoso e rejeito esteril) apresentam urn alto grau 
de confusao entre si e com a classe area urbana, quando classificadas 
automaticamente (Kux e Valeriano, !981). 

Por esse motivo o algoritmo SELECAO DE ATRIBUTOS foi uti_ 
lizado na escolha dos quatro atributos que melhor discriminam estas 
3 classes de uso da terra. Para isto, amostras de treinamento destas 
3 classes foram extraTdas no vTdeo do IMAGE-100 e fornecidas ao alg£ 
ritmo SELECAO DE ATRIBUTOS. 


*Distdnoia JM - distancia estatistica que , vaviando de zero a 2, refle_ 
te o grau de dissimilar idade entre duas distribuiqoes , sendo zero o va_ 
lor enoontrado entre distribuiqoes totalmente semelhantes , e o valor 
2 para distribuiqoes sem "overlapp" (Swain and Davis , 1978, p. 170- 
174). 


2.4 - CLASSIFICACAO DOS DADOS MSS-LANDSAT 


, Classificar uma imagem atraves de analise digital consi£ 

te em atribuir os "pixels" da cena a determinadas classes dentro de uni 

criterio estatTstico ou deterministic. A abordagem util izada neste tra 

balho foi a classificacao determinTstica nao-supervisionada (real izada 

pelo proprio analisador de imagens), cujo resultado orienta a localize 

> 

cao de amostras de treinamento para uma classificacao estatTstica s_u 
pervisionada (classificador hTbrido, segundo Swain and Davis (1978). 

2.5 - CLASSIFICACAO NAO-SUPERVISIONADA 


Para efetuar esta classificacao foi utilizado o algorit 
mo K-MEDIAS. Este algoritmo discrimina nuvens de adensamento de popula. 
Coes dentro do espaco de atributos utilizado. Os centros destas nuvens 
sao local izados atraves de urn procedimento de medias errantes, e as de 
;nais cel as do espaco de atributos sao atribuTdas as classes dentro de 
urn criterio de menor distancia euclidiana ao centro da nuvem. 

0 resultado deste procedimento e a estratifi cacao da cei 
na em classes que agrupam "pixels" de padroes espectrais semelhantes. 
Este resultado e util na selecao de amostras de treinamento numa cla£ 
sificacao estatTstica, uma vez que ele diminui a possibil idade de for 
necer ao treinamento classes com comportamento espectral bimodal (como 
o algoritmo utilizado na classificacao supervisionada supoe que a dis_ 
tribuicao do comportamento espectral da classe e gaussiana, tal erro 
compromete todo o procedimento). 

Oito classes foram obtidas da analise nao-supervisionada. 
Este resultado foi comparado a verdade terrestre, e amostras de treina 
mento foram selecionadas em areas com cobertura do solo homogenea e 
classificacao nao-supervisionada uniforme. 


Os atributos selecionados em ambos os criterios (maior 
distancia J.M. media e maior distancia J.M. minima) foram: 

MSS-5. 

» 

MSS-6. 

RAT 7/4. 

RAT 6/5. 

♦ 

* As distancias J.M. entre as classes e os respectivos li 
mites inferior e superior da probabil idade de erro de classificagao e^ 
contram-se na Tabela 1. 


TABELA 1 

CALCULO DAS DISTANCIAS J-M 


'''""\CANAIS 

CLASSES"\^ 

MSS 5, MSS 6, RAT 7/4, RAT 6/5 ■ 

1 e 2 

DIST.J-M = 1,1508 
LIM. INF. = 4,73050E-02 
LIM. SUP. * 0,21229. 

1 e 3 

DIST.J-M = 1 ,9572 
LIM. INF. = 1 ,1 4381 E-04 
LIM. SUP. = 1 ,0694£-02 

2 e 3 

DIST.J-M = 1,3342 
LIM. INF. = 2,85217E-02 
LIM. SUP. = 0,16646 

1, 2 e 3 

DIST.J-M MEDIA = 1,4807 
LIM. SUP. TOTAL = 0,38944 

DIST.J-M MINIMA = 1.1508 (entre classes 1 e 2) 
LIM. INF. = 4.73050E-02 
LIM. SUP. = 0,21229 


CLASSES: 1. Rejeito piritoso 

2. Rejeito esteril 

3. Area urbana 



2.6 ~ CLASSIFICACfiO SUPERVISIONADA 


Apos a selegao das areas para a amostragem de treinamen. 
to, foram caracterizados os vetores medios e as matrizes de covariancia 
nos quatro canais de doze classes de uso da terra, abaixo discriminadas: 

- Duas classes de depositos de rejeito piritoso, uma corresponded 
do a depositos recentes pouco oxidados e a outra, a depositos an 
tigos relativamente mais claros. 

- Duas classes de depositos de rejeito esteril, as quais difereji 
ciam rejeitos produzidos por maquinas diferentes e, consequent^ 
mente, tamanhos de cones de deposigao diferentes. As areas coin 
grandes cones (associados a oriel? tagao da lavra em relagao ao 
azimute do sol) apresentam-se mais escuras na iniagem devido ao 
sombreamento gerado pelos proprios cones. 

- Tres classes de areas agrTcolas: uma referente a pastagens e 
areas cultivadas com cobertura vegetal na epoca da passagem do 
satelite, e duas classes correspondentes a areas preparadas p^ 
ra plantio, com solos apresentando diferengas espectrais detecta 
das na classificagao nao-supervisionada. 

- Duas classes de areas florestais: uma abrangendo matas naturais, 
eucaliptais desenvol vidos e capoeiras, e outras caracterizando 
apenas eucaliptais novos ou rebrotas tambem em desenvolvimento. 

- Tres classes de areas urbanas: uma correspondcnte a cidade pr£ 
priamente dita, com tonalidade cinza-medio e texturizada, outra 
a zona periferica da mancha urbana, correspondente a areas de 
baixa densidade habitacional , e a terceira caracterizando as in 
dustrias. 

Os parametros estatTsticos de cada classe foram fornec£ 
dos a urn classificador bayesiano de imagens digitais (algoritmo MAXVER) 


is-rr 


nr 
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que atribui cada "pix&V da cena a uma das classes dentro de um critS 
rio de maxima verossimilhanga. Na representagao tematica deste resulta 
do, as classes afins sao apresentadas em um unico terna, gerando assim 
um mapa tematico que apresenta as seguintes classes: 

- rejeito piritoso, 

- rejeito esteril, 

- area agrTcola, 

- area florestal, 

- area urbana. 

Apresenta-se o resultado deste procedimento na Figura 1, 
uma compos icao colorida (MSS 5 em cyan e MSS 7 em vermelho) da mesrna 
area na Figura 2 para comparagao e o produto final na Figura 3. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 



Fig. 1 - Classificagao automatica do uso do solo 
(Modulo Criciuma - passagem 24/04/78). 

Legenda: Azul escuro - Rejeito piritoso 
Amarelo claro - Rejeito esteril 
Amarelo escuro - Area agrTcda 
Verde - Area florestal 
Lilas - Area urbana. 
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Fig. 2 - Composigao colorida da area de estudo. 

MSS 5 - cyan 
MSS 6 - vern.elho 

(aos dois canais foram aplicados 
"scaling" x 2) 



A seguir Figura. 3. 
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3. ESTIMATIVA DA EXATIDAO DE CLASSIFICACAO 


Exatidao de classificagao de nrnpas tematicos e uma medj. 
da do quanto o mapa representa a realidade, ou seja, o percentual de 
pontos corretamente classificados presentes num mapa tematico. 

Para efetuar esta estimativa, realizou-se a interpret^ 
gao de fotografias aereas da area de estudo, na eseala 1:45.000, cujo 
sobrevoo foi realizado em fins de agosto de 1978, A area de estudo e 
coberta por 21 fotografias arranjadas em tres faixas com recobrimento 
lateral de 10 a 40% e longitudinal de 60%. Supos-se que o mapa deriva 
do das fotografias aereas foi 100% correto e denominado mapa de verda_ 
de terrestre (MVT). 

Uma amostragem de pontos foi selecionada ao acaso no MVT. 
Estes pontos foram local izados no mapa derivado de dados LANDSAT, e o 
resultado da comparagao entre as duas classificagoes foi computado. 

3.1 - ELAB0RACA0 DO MVT 


Como a final idade un i ca desta fase do trabalho e obter 
dados de verdade terrestre para comparagao com o resultado da anal ice 
digital dos dados MSS-LANDSAT, as classes discriminadas neste mapa apr^ 
sentam urn grau de detalhamento proximo ao das classes obtidas dos d^ 
dos LANDSAT. 


• Foram extraTdos em "overlays", juntamente aos principal’s 
rios e estradas, os 1 iniites das seguintes classes: 

1) rejeito piritoso, 

2) rejeito esteril, 

3) area agrTcola, 

4) area florestal , 

5) area urbana, 
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6) area alagada, 

7) agua, 

8) solo exposto. 

Em seguida, estes "overlays" foram colocados em forma de 
mosaicos utilizando como base cartografica 12 pontos de controle (4 por 
faixa), extraTdos da carta topografica CRIClOMA elaborada pelo IBGE na , 
escala 1:50.000 e‘ corrigidos para a escala das fotografias. Dada a baj_ 
xa densidade de pontos de controle, o mosaico gerado nao e controlado 
cartograficamente. 


Apresenta-se o resultado deste procedimento na Figura 4. 
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3.2 - AMOSTRAGEM PARA ESTIMATIVA DA EXATIDAO 


Uma vez que as classes 6, 7 e 8 do MVT tern baixTssinia ex 
pressao em area, estas nao foram consideradas neste procedimento. 

Cinquenta pontos foram escolhidos ao acaso no MVT sobre 
cada uma das seguintes classes: 

- rejeito piritoso, 


- area agrTcola, 

- area f 1 oresta.1 , 

- area urbana, 

e trinta pontos sobre a classe rejeito esteril (ja que sua area tambem 
e pequena). 


A local izacao destes pontos foi assinalada numa tela s£ 
bre a qua! havia sido projetado um diapositivo referente a Figura 2, 
ampliado a uma escala proxima a do MVT. As fotografias aereas foram 
utilizadas para auxiliar esta transposicao. 


Em seguida, foi projetado o diapositivo da Figura 1, p£ 
sicionado no mesmo local da projecao anterior (com uso de pontos de co£ 
trole visualmente identificaveis nos dois diapositivos) , e a classify 
cacao de cada ponto foi observada gerando uma matriz de confusao, apre 
sentada na Tabela 2. 
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4. CONCLUSflO 


Um teste de hipoteses num nTvel de significancia de 0,95 
perirri te afirmar que as duas primeiras classes (area florestal e area 
urbana) estao discriminadas com uma exatidao superior a 85# (Ginevan, 
1979). A classe area urbana, apesar de incluir areas dos dois tipos de 
rejeito, apresenta uma pequena margem de omissao, principalmente na pe 
riferia das cidades, confundindo-se com a classe area agrTcola . As duas 
classes de rejeito nao foram classif icadas de modo satisfatorio, segun 
do o teste realizado. 

Mesmo considerando que o risco de cometer um erro de do 
cisao do tipo II (este rejeitou um mapa com exatidao superior a 0,85, 
no caso um mapa de exatidao = 0,90) para este tamanho de amostra e de 
75%, nao se pode recomendar processamento automatico de dados MSS- 
LANDSAT para mapear as classes de rejeito de carvao. 

Recomenda-se para a execucao da parte final do projeto, o 
uso de tecnicas de interpretaqao visual apl icadas a composiqoes colon 
das para a extracao das areas ocupadas por rejeitos de carvao. Estas 
tecnicas aliadas a realcamento de imagens e informacoes auxiliares pa_ 
recem bastante viaveis, uma vez que estas classes ocupam areas restn 
tas e bem diferenciaveis visualmente. 

As demais classes deverao ser identif icadas automatic^ 
mente, apenas com a precaucao de extrair visualmente os limites das 
manchas urbanas. 

Finalmente, adianta-se que, segundo observaqoes em modu^ 
los costeiros (Ararangua, Laguna e Imariu (Valeriano e Pereira, 1983)), 
dever-se-a utilizar tecnicas de analise visual de dados MSS-LANDSAT na 
separagao entre as classes areas florestais e areas alagadaSa 
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